arXiv:1506.08229v2 [nlin.SI] 24Jul2015 


Dimer with gain and loss: Integrability and "PT-symmetry restoration 

I. V. Barashenkov 

Centre for Theoretical and Mathematical Physics, University of Cape Town, 

Rondebosch 7701, South Africa and Joint Institute for Nuclear Research, Dubna, Russia 

D. E. Pelinovsky 

Department of Mathematics, McMaster University, Hamilton ON, 

Canada, L8S fKl and Department of Applied Mathematics, 

Nizhny Novgorod State Technical University, Nizhny Novgorod, Russia 

P. Dubard 

Department of Mathematics, University of Cape Town, Rondebosch 7701, South Africa 

A "PT-symmetric nonlinear Schrodinger dimer is a two-site discrete nonlinear Schrodinger equation 
with one site losing and the other one gaining energy at the same rate. In this paper, two four- 
parameter families of cubic "PT-symmetric dimers are constructed as gain-loss extensions of their 
conservative, Hamiltonian, counterparts. We prove that all these damped-driven equations define 
completely integrable Hamiltonian systems. The second aim of our study is to identify nonlinearities 
that give rise to the spontaneous PT-symmetry restoration. When the symmetry of the underlying 
linear dimer is broken and an unstable small perturbation starts to grow, the nonlinear coupling of 
the required type diverts progressively large amounts of energy from the gaining to the losing site. 
As a result, the exponential growth is saturated and all trajectories remain trapped in a finite part 
of the phase space regardless of the value of the gain-loss coefficient. 
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I. INTRODUCTION 

The nonlinear Schrodinger dimer is a code name for the discrete nonlinear Schrddinger equation defined on a lattice 
consisting just of two sites: 


ill + v = F{u, u*, V, V*), 
iv + u = G{u, u*, V, V*). 


Here and in what follows, the overdot denotes the derivative with respect to t. 

The dimer is one of the simplest (and hence most heavily used) models of a coupled bimodal structure. In optics, 
the system 0 describes the directional coupler — a pair of parallel waveguides coupled through their evanescent fields. 
In this context, u and v are the complex amplitudes of stationary light beams in the waveguides and t measures the 
distance along their parallel cores [TH3!. 

When the same system is employed in the studies of the Bose-Einstein condensates, u and v stand for the amplitudes 
of the mean-field condensate wave functions localized in the left and right well of a double-well potential [J] (or of 
their symmetric and antisymmetric combinations 0 )- In this application, t has the meaning of time. 

The nonlinear Schrodinger dimers were also utilised in the solid state physics [SHE] and in the context of electric 
lattices [5]. 

Typically, the dimer arises as an amplitude equation; that is, u and v represent slowly changing amplitudes of some 
oscillatory variables x and y: 

x{t) = -b ..., 2 /(r) = + z;*(t)e-“^ -b .... 

Here t = e^r, is a small parameter, and the dots stand for small anharmonic corrections. The invariance of the 
optical or atomic structure with respect to translations in r is inherited by the amplitude equations as the invariance 
under simultaneous phase shifts in u and v. In other words, physically meaningful nonlinearities have to satisfy 

F{e’^*u, e-^^u*,e^^v, =F’^F{u, u*, v, v*), 

G{F^u,e-'^u*,e^'^v,e-^'^v*) =e^^G{u,u* ,v,v*) 


for any real cj). 

Another property of the dimer dictated by physics, is conservation of energy which gives rise to an underlying 
Hamiltonian structure. There are two main types of Hamiltonian formulations admitted by the dimers. One class of 
Hamiltonian dimers has the straight-gradient form 

.du 8% .dv 8% 

^~di ^ 8^' ^di ^ ^ ’ 

Here % is the Hamilton function while the canonical coordinate-momentum pairs are (m, u*) and (u, v*). Alternatively, 
u can be paired with v*, while u* play the role of momentum conjugate to v. 

.du 8^ . dv 8S) 

dt dv* ’ * dt 8u* 

(This time we use a different notation for the Hamilton function to keep the treatment of the two cases separate.) In 
what follows, we are referring to equations Q as the cross-gradient systems. 

The last requirement is that of the left-right symmetry (the parity symmetry) of the system. This requirement 
arises if the two elements making up the dimer (two waveguides or two potential wells) are identical. Mathematically, 
it reduces to the invariance under the permutation of u and u in Q: 

F{u, u*,v, V*) = G{v, V* ,u, u*). (5) 


If the two elements are not identical — for example, if one channel is dissipative while the other one draws energy 
from outside — the system may be still invariant under a relaxed form of the left-right permutation known as the 
parity-time (VF) symmetry [3]. Mathematically, the discrete Schrodinger equation is said to be T^T-symmetric if it 
is invariant under the product of V and T transformations. Here the V operator swops the two elements around. 


r 




while the T represents the effect of time inversion on the complex amplitudes: Tu{t) = u*{—t), Tv{t) = v*{—t). 
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The current upsurge of interest in the T^T-symmetric systems is due to the fact that they can strike the balance 
between the gain of energy in one channel and loss in the other. In the optical context, the T^T-symmetric dimer 
describes a waveguide with power loss coupled to a waveguide experiencing optical gain miiiinHis]- In the matter- 
wave setting, the T^T-symmetric system is formed by two quantum states, with one state leaking and the other one 
being fed with particles [Hlig. 

Since the T^T-symmetric systems have channels for the energy exchange with the environment, they are commonly 
thought to occupy a niche between dissipative and conservative systems. It was therefore met with surprise when 
some linear M and nonlinear m T^T-symmetric systems were found to possess Hamiltonian structure. In particular, 
there are Hamiltonian T^T-symmetric dimers; one example was produced in HZl: 

iii + v + (\u\'^+ 2\v\'^)u + v^u* = iju, 
ii) + U+ (|up -I- 2 |mP)u -I- u^v* = —i^v. 

The two terms in the right-hand sides of (§ account for the gain and loss of energy, with 7 > 0 being the gain-loss 
coefficient. 

Another Hamiltonian T^T-symmetric system was identified in |18j : 


iu + v + \u\'^u= i^u, 
iv + u + \v\'^v = —i"fv. 


( 7 ) 


Because of its ubiquity in physics nnniHiiis], equation Q is occasionally referred to as the standard dimer. 
Finally, the Hamiltonian model 


iu + V — {2ai + a 2 )\u\'^v — 2ai\v\'^v = i^u, 
ii) + u — (2ai -\- a 2 )\v\'^u — 2ai\u'\^u = —i^v, 


( 8 ) 


was discovered outside the domain of the VT -symmetry — as a by-product in the search of integrable equations [ 20 ] . 
Here ai and 02 are arbitrary real coefficients. 

The availability of the Hamiltonian structure is a fundamental property of a dynamical system. This property by 
itself implies the conservation of phase volume and hence some degree of regularity of motion. But in the presence of 
additional first integrals it allows to establish an even higher level of regularity, namely, the Liouville integrability. 

The first aim of this paper is to show that any cubic T^T-symmetric phase-invariant dimer obtained as a VT- 
symmetric extension of the conservative dimer ([^ or Q, is a Hamiltonian system. By determining an additional 
integral of motion independent of the Hamiltonian, we establish the complete integrability of all these systems. 

Another topic pursued in the present study concerns the phenomenon of T^T-symmetry breaking — one of the 
experimentally accessible properties of physical systems with gain and loss |12[ I21I123| . The spontaneous symmetry 
breaking occurs in linear T^T-symmetric systems as the gain-loss coefficient is increased beyond a critical value 
7c. This exceptional point separates the symmetric phase (7 < 7c), where all perturbation frequencies are real, 
and the symmetry-broken phase (7 > 7c), where some frequencies are complex and the corresponding modes grow 
exponentially. 

When the input power in the physical structure is low — or, equivalently, when the initial conditions of the 
corresponding mathematical model are small — all nonlinear effects are negligible and the system follows the linear 
laws. In particular, small initial conditions in the symmetry-broken phase trigger an exponential growth. However 
as the resulting solution reaches finite amplitude, the nonlinear coupling terms kick in. These terms may channel the 
power from the site where it is gained, to the site where it is lost. The higher is the power gained, the larger portion 
of it is channeled to the disposal site by the nonlinear coupling. In systems where this mechanism is at work, the 
exponential growth is arrested and all escaping trajectories are sent back to the finite part of the phase space. The 
VT symmetry becomes spontaneously restored. 

The classihcation of integrable T^T-symmetric dimers with the nonlinearly restored 7^T-symmetry, is the second 
objective of our study. 

The outline of the paper is as follows. 

In section |ll] we present a four-parameter family of the VT -symmetric cross-gradient dimers; all these systems are 
Hamiltonian in their original u and v variables. In a similar way, a four-parameter T^T-symmetric extension of the 
straight-gradient dimer is introduced in section [ml 

The straight-gradient systems do not admit the Hamiltonian formulation in terms of u and v — except when the 
straight-gradient system is cross-gradie nt at the same time, or when it is gauge-equivalent to a cross-grad ient system. 
(We identify such dual cases in section HI ) Nevertheless, transforming to the Stokes variables (section IVI we can 


describe all trajectories of the straight-gradient dimer and elucidate the geometry of its phase space. 
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In the subsequent three sections we determine the canonical coordinates for the general straight-gradient dimer 
and reformulate it as a Hamiltonian system. Three complementary subfamilies of the straight-gradient models are 
considered (sections [V VI and |VII[ ) ■ 

Section [VlII| is concerned with the trajectory confinement and 'PT-symmetry restoration. We identify broad classes 
of nonlinearities capable of suppressing the exponential blowup regimes — both within the cross-gradient and straight- 
gradient families. 

Section |IX| offers examples of simple oscillatory systems with the amplitude equations in the form of cross- and 
straight-gradient dimers. 

Finally, in section!^ we summarise mathematical results of this study and discuss their physical implications. 


II. CROSS-GRADIENT PT-SYMMETRIC DIMER 


A general cross-gradient dimer (|3, complying with the phase invariance Q, with no gain or loss, with cubic 


nonlinearity, permutation property and linear part of the form Q, is defined by the Hamiltonian 


^0 = -{\u? + kP) + W{u,v). 


(9) 


Here IF is a U(l)-invariant real quartic polynomial in u, v and their complex-conjugates, which is symmetric with 
respect to the u ^ v permutations. The most general quartic polynomial with these properties can be written as 

W = ai(|up -I- \v\'^)'^ + a 2 |wnup -I- a^{u*v + uv*){\u\'^ + jup) -I- a 4 {u*v + uv*Y , (10) 


where 01,021 and 04 are real coefficients. (See Section 2 in [H]). 

For any set of o’s, this gainless lossless system admits a straightforward PT-symmetric extension 

.du dSj . dv dSj 

* df “ 

where the Hamilton function Sj is different from i^o in just one term: 


( 11 ) 


Sj = —(|mP -I- |up) -|- W{u,v) + ij{uv* — u*v). 


( 12 ) 


Here 7 > 0 is the gain-loss coefficient. Substituting the expression with W 
a four-parameter family of Hamiltonian PT-symmetric cubic dimers: 


as in ( 10 ) in 


equations 0 , we obtain 


iii + v — i^u = oadup -I- 2\vf')u -\- a^v'^u* + 2a4U^v* + [(2oi -|- 02 -|- 2a4)|up -I- 2oi|wp] v, 

iv + u + i^v = 03 ( 2 |up -I- lup)?; -I- a^u^v* + 2a4v‘^u* + [(2oi -1-02-1- 2o4)|up -I- 2 oi|mP] it. 


(13) 


A particular case of (13) is the system §. This is selected by letting 03 = — 1 and oi = 02 = 04 = 0. Another 
special case is the dimer Q; this model corresponds to 03 = 04 = 0. The Hamiltonian structure of these two 
particular systems has been determined earlier [171 120|. 


Despite the seeming complexity, all trajectories of (13) admit a simple analytic description. We define the Stokes 
vector TZ — \X +}Y + kZ, where 


X = u*v + v*u, Y = i{u*v — v*u), Z = |iip — |i;p. 


(14) 


Note that the length TZ 


\JX'^ -|- -|- Z^ of the Stokes vector has a simple expression in terms of u and v: 


n = \u? + \v?. 


Transforming to X, Y, and Z, equations (13) give a dynamical system in three dimensions: 


V = 0, 

Y = -2Z + 2a3XZ + 4oiZP, 
z = 2Y + 2-fn - 203 vr - (02 -f 4oi)rp. 


(15a) 

(15b) 

(15c) 


Equation (15a) implies that for any selection of Oi, 02 , 03 , and 04 , the cross-gradient dimer (131 has two independent 


integrals of motion: X and Sj. Accordingly, the Hamiltonian system (13) is Liouville-integrable 
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All trajectories lie in parallel vertical planes X = Xq, where Xq is an arbitrary constant. The form of the trajectories 
is determined by setting X = Xq in the equation Z) = const, with Sj being the Hamiltonian (12) expressed 

in the Stokes variables: 


Y - 7 ) + ^ {aiTZ + asXo - 1) = C. 


(16) 


Here C is a constant of integration, 
trajectories on the X = Xq plane. 


Equation (16) with TZ 


y/Xl +Y^ + Z^ defines a one-parameter family of 


III. STRAIGHT-GRADIENT PT-SYMMETRIC DIMER 


The conservative straight-gradient dimer with general cubic nonlinearity, left-right symmetry and phase invariance, 
has the form 


where 


, du dV. . dv dH 


(17) 


Ti. = —{uv* + u*v) + W{u, v), 


(18) 


and W is the four-parameter quartic polynomial: 


W = /3i(|up -I- \v\'^Y + P 2 \u\^\v\'^ + I33{u*v + uv*){\u\‘^ -I- |up) -I- j3i{u*v + uv*Y. 


(19) 


Here /3i, /32, /^s, and /34 are real coefficients. This is the same quartic as in the previous section; we have just switched 
from the a- to the /3-notation to emphasise that we consider a totally new family of models. 

The system ([T7|)-([l^ with /3i = — i, /32 = 1, and P 3 = = 0 is known as the N = 2 discrete self-trapping equation 

[5]. When /33 = /34 = 0 while /3i = —p/2, /32 = p — 1 with p a real coefficient, these equations constitute the spatially 
homogeneous version of the Aceves-Wabnitz coupled mode system for the nonlinear optical grating |25j . A particular 
case of this (/3i = 0, /32 = —1) is the spatially-independent Thirring model (a theory of self-interacting spinor field) 
[55H1H1- Another case related to spinors is /34 = — j, /3i = /32 = /da = 0; this system derives from the one-component 
Gross-Neveu model [23|29]. The system with /32 = 2, /34 = — |, /3i = /33 = 0 is related to the spinor theory with the 
pseudoscalar interaction [2H1130] • 

The T^T-symmetric extension of the general straight-gradient dimer (|17[) has the form 


du 

i— - I'yu = 

dt 


cm 

c/u*’ 


dv 

l— -h I'^V 

dt 


dU 

dv*' 


Evaluating the partial derivatives using (18) and (19), these equations become 


iii + v- i'yu = [2/3i|up -I- (2/3i Y P 2 + 2/34)|up] u + 2j3iv'^u* + jd^u^v* + jd^ (2|up -I- |up) v, 

iv + U + i'yv = [2,01 |up -|- (2/3i -I- ,02 + 2 / 34 )|up] v + 2fdiU^v* + Pqv'^u* -|- Pq (2|up -I- |up) u. 


( 20 ) 


( 21 ) 


The standard dimer ([^ is a special case of (21). This is selected by taking fdi = — ^2 = 1, and jd^ = jd^ = 0 in 

equations (21). There is an extensive literature on mathematical aspects of this model [H [H [3IHS7 ■ 

The couples {u,u*) and {v,v*) do not form pairs of canonically conjugate variables. That is, the straight-gradient 
dimer (21) does not admit a Hamiltonian formulation in terms of the original complex coordinates — except when 
the straight-gradient dimer is cross-gradient at the same time. The necessary and sufficient condition for the equation 
0 to have a representation 0 with some ij, is 


d 

( ■ 

d 1 


du* 

('^“+a?) 

1 * 

II 



For H of the form (18), this condition translates into 


du*"^ dv*^' 


Substituting the quartic polynomial (19) for W gives, finally, /3i = P 4 . 
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The choice /3i = 13^ ensures the existence of a (complex) function io such that 


i^u ■ 


-^7^! ■ 


dn 

dS) 

, , m 

5Jo* 

( 22 ) 

du* 

~ dv*' 

+ a? 

dv ’ 

dn 

dSj 

, , dH 

9J5* 

(23) 

^ ' 

~ 

n. +^ = 

du 


The necessary and sufficient condition for the function Jo in (22) to be real, is given by 

d 


dv 


dn 


i^u - 


d 


du* I dv* ( du 


dn 


In a similar way, the necessary and sufficient condition for Jo = Jo* in equation (23) is 

d 


dn 


d 




dU 


For n of the form (18), each of these two conditions amounts to 

d^W d^W 


dv*du du*dv^ 

which gives 2/3i (32 = 2 / 34 . Using the previously established condition [3i = j3i, this relation reduces simply to 

/?2 = 0 . 


Thus, the straight-gradient dimer (21) with /3i = Pi and /32 = 0 is, at the same time, a Hamiltonian system with the 
cross-gradient canonical structure (11). The corresponding Hamilton function Jo is determined by simple integration: 


^ = -(|up + lup) + + 2Pi{\u\^ + 




vf)(uv* + u*v) + —{uv* 


+ vu*Y ni'^iuv* — vu*). (24) 


The Hamiltonian (p4|) is of the form ([T^, ((Tol) with = 04 = 02 = 0, and 0:3 = 2 / 34 . 


In fact, the class of straight-gradient dimers admitting the Hamiltonian formulation is even wider. Let uq and vq be 
a solution to the straight-gradient equations (211 with P 2 = 0 and generic Pi, Pi (that is, Pi not necessarily coinciding 
with Pi). The gauge transformation 


Uq = e^^u, Vq = e^'^v 


with 


^{t)=2{Pi-pi) f{\u\^ + \v\^)dT 
Jo 

generates functions u and v which satisfy the straight-gradient dimer equations with P 2 = 0 and Pi set equal to 
iii + v — i^u = 2 Pi (|up -|- 2 |up) u + 2 PiV^u* + Psu^v* + P 3 ( 2 |up -|- |up) v, 
iii + U + i"fv = 2Pi (|vp -|- 2 |mP) v + 2PiVpv* + Pov'^u* + P 3 (2|up 


lup) u. 


(25) 


As we already know, this system has a cross-gradient Hamiltonian formulation with the Hamilton function ( |24[ ). 
Therefore, the straight-gradient dimer with P 2 = 0 and any value of (Pi — Pi) is gauge-equivalent to the cross-gradient 
Hamiltonian system where the Hamilton function Sj is as in ([24|). 


In what follows, we uncover the Hamiltonian formulation of the straight-gradient system (21) with P 2 ^ 0. This 
class will include, in particular, the standard dimer 0 (which has P 2 = 1 ). 


IV. PHASE SPACE OF STRAIGHT-GRADIENT DIMER 


The generic {P 2 ^ 0) straight-gradient dimer does not admit the Hamiltonian formulation in terms of the original 
u and V variables. In order to determine the canonical pairs of coordinates, we transform it to the Stokes variables 
(14). Equations @ give a dynamical system in three dimensions: 


X = -P 2 YZ, 
Y = -2Z + {P 2 + 4pi)XZ + 2P3ZTZ, 
Z = 2Y + 277^ - 4PiXY - 2/337^U. 


(26a) 

(26b) 

(26c) 
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We also note an equation for the length of the Stokes vector that follows from the system (26): 

h = 27 Z. 


(27) 


Replacing t with 77. as a new independent variable, and using d/dt = TZd/dTZ, equations (26a) and (26b) become a 
linear nonhomogeneous system 




dY 


27^ - {132 + 4/34)^ = -2 + 2/3377. 


(28) 


Our strategy will be to determine the general solution of ( |28[ ) including two constants of integration. These “constants 
of motion” of the system (28) will serve as the two first integrals of the original three-dimensional dynamical system 
(26). To obtain the Hamiltonian formulation of the dimer (21), one of these will be appointed as the Hamiltonian 
and the other one as a canonical coordinate. 

As in the case of the cross-gradient dimer, the existence of two independent conserved quantities along with the 
availability of the Hamiltonian structure will imply that the straight-gradient dimer (21) is Liouville integrable. 

It is convenient to introduce the quantity 

w2 = /32(/32 + 4/34). (29) 

The form of the solution of the system (28) depends on whether is positive, negative or zero. We consider these 
three cases separately. 


Assuming that > 0, the general solution of (28) is 


X = A cos I —77 

.27 


i? sin I ^TZ 


27 


P 2 


27 


A sin —77 I — B cos —77 


27 


2 - 2/3377 
/32 + 4/34 ’ 
47^3 

O 5 


(30a) 

(30b) 


where A and B are constants of integration. Treating 77 as a parameter, and supplementing (30al-(30b) with the 
formula 


Z = ±\/372 - X2(77) - Y^iTZ), 


(30c) 


these equations provide explicit parametric expressions for trajectories of the dimer: X = X{TZ), Y = Y{TZ), 
Z = Z(77). 


Denoting p = VA^ -I- 73^, equations (30a) and (30b) give 

2 / 3377 - 2 ^ ^ 

/32 + 4^4 


P=\X + 


-I (T 


_ 47/33 \ 


2 

2 j ■ 


(31) 


Equation (31) with 77 = y/X"^ +Y'^ -\- Z"^ is an implicit equation of the surface on which all trajectories lie. The shape 

(32) 


of this surface in the (AT, Y, Z) phase space depends on the value of the parameter 

2/33 


/32 -I- 4/34 

Figs[^(a), (b), and (c) depict the surface with \a\ < 1, |cr| > 1, and \a\ = 1. 


In the case where uP < 0, we define > 0. The general solution of (28) is then 


X = C exp ( 7^77 ) -I- 77 exp ( —;^77 


V27 


27 


^ P 2 


27 


C exp 7 ;—77 I — D exp — —TZ 


27 


2 - 2/3377 

/32 + 4/34 ’ 

47/33 

,/2 ’ 


where C and D are two constants of integration. 

Expressing 4C77 in terms of X, y,77 and denoting it siS 2 P^, where Si = signC and S 2 = sign 77, we obtain 

47 ^ 3 \^ 


S 1 S 2 P = X -I- (t77 - 


02 + 4/34 






(33a) 

(33b) 

(34) 


Here p> 0 and S 1 S 2 is either 1 or —1, depending on whether CD > 0 or CD < 0. The solution surface (34) made up 
by all trajectories of the dimer, is shown in Fig[^ (d), (e), and (f) for \a\ < 1 , \a\ > 1 , and \a\ = 1 , respectively. 
Finally, the case = 0 is considered in section [VIl] below. 





































FIG. 1. The surface (311 with > 0 (a, c, e) and < 0 (b, d, f). The parameter |(t| is smaller than 1 in (a,b); greater than 
1 in (c,d), and equal to 1 in (e,f). The value of /32 is 1 in (a,c,e) and -1 in (b,d,f). The value of /Is is 2 in (a), 1 in (b), 5 in (c), 
9/2 in (d,e), and 7/2 in (f). In all panels, 7 = 1 and = 2. The surface parameter = 7 in (a,c,e) and SiS 2 p^ = in 

(b,d,f). We also show a few trajectories of the system (261 lying on each surface. 


V. HAMILTONIAN STRUCTURE OF STRAIGHT-GRADIENT DIMER: > 0 


It is convenient to write the explicit solution (30a)-(30b I in the complex form 


W P2 + 4/54 


( 35 ) 
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where 


pe * = [A —iB)ex.p\i—TZ 


OJ 


(36) 


and p = y /was introduced in the previous section. Assume p > 0; in this case, equation (36) serves as the 
definition of 6. (We will examine the possibility p = 0 in subsection VD below.) Using equation (35), 0 is expressible 
as a function of X, Y, and TZ: 


0 = arctan 


oj X (jTZ — ‘2/32^ 
(32 Y - 47/33W 


-2 


-2 


(37) 


The definition (36) implies that the complex quantity 


is a constant. Therefore its argument 


pexp — (—270 — loTZ) ^ = A — iB 


H = -2j0 - w7^ 


(38) 


is a first integral of the system (26). We shall demonstrate that H can serve as a Hamilton function for the three- 
dimensional system (26) and the underlying nonlinear Schrodinger dimer (21). 

One more integral of motion, p, is given by (31). Choosing p and 0 as a pair of coordinates in the phase space, we 
require that TZ (and hence H) be a function of p, 0, and Pg — but not depend on Pp [THl. Irrespectively of how we 
define Pp, the Hamilton equation 


P = 


dH 


(39) 


will then reproduce the equation for p: p = 0. Another consequence of requiring dH/dPp = 0, is that the variable Pp 

(40) 


will not participate in the dynamics and the Hamilton equation 

dH 


'dp 

will decouple from the rest of the system. 

At this point, we note that, since 0 is only different from —^TZ by a constant, equation (27) gives 

0 = —LuZ. 

We should choose the momentum Pg (denoted P for brevity) in such a way that the canonical equation 

dH 


(41) 


0 = 


dV 


(42) 


reproduces equation (41). For the Hamilton function of the form (38), equations (42) and (41) are equivalent if 

dTZ 


dP 


= Z. 


(43) 


Once the equation (42) is satisfied, the conjugate equation 

dH 


V = - 


de 


(44) 


will be verified automatically. Indeed, since H does not depend on Pp while p = 0, there are only two nonvanishing 
terms in the derivative H: 

■ dH- dH- 


Substituting for 0 from (42), the relation H = 0 implies Eq.(44). 
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Thus all we need to do in order to put the canonical structure in place, is to identify the canonical momentum 


V = Pg ensuring the validity of equation (43). Before proceeding to the unveiling oi Pg, a technical remark is in order. 


Equation (35) implies that the coordinate Y can be expressed solely in terms of p and 9 (rather than p, 0, and TZ)\ 


47/33 w 

Y = —^ + —pcosO. 
d 2 


(45) 


Therefore 9 may be thought of as the spherical polar angle in the frame of reference where Y is the vertical coordinate. 
In what follows, we also introduce (an analogue of) the azimuthal angle on the (X, Z)-plane. This construction will 
be based on the following decomposition of the coordinate X stemming from (35): 

(46) 


X = X — aTZ. 

Here we have isolated a term that is expressible entirely in terms of p and 9: 

2 


X = 


h + 4/34 


psind. 


(47) 


Our construction of the canonical momentum Pg depends on whether the parameter a defined in (32) is smaller, 
greater, or equal to 1 in magnitude. 


A. Hyperbolic case: |(j| < 1 


Assuming |cr| < 1, we dehne 


= 1 - > 0 . 


The corresponding subfamily of models includes, in particular, the standard dimer ([^ — for which = 1 and tr = 0. 
(We note that the Hamiltonian structure of the standard dimer was elucidated earlier [H].) 

Inserting the decomposition (46) in the identity X^ + Y'^ + Z'^ = P? yields 




(48) 


where we have substituted for Z from (43) and defined 


A2 


(49) 


Since x and r are expressible in terms of p and 9 only (i.e., are independent of P), equation (48) is an ordinary 
differential equation for P{V). To define 7^, it is sufficient to pick one solution of this separable equation. We choose 

P{V) = ^ cosh(7 


the formula (43) gives then 


so that 


Z = r sinh(A7^), 


Pg = V = — arcsinh ( — 


cr 

- 

(50) 

), 

(51) 

f z\ 

(52) 

7 ■ 


Together with the definition of the coordinate p in (31), the coordinate 9 in (37), the definition (52) completes the 
set of three canonical variables. The fourth coordinate, Pp, can be reconstructed from equation (40) by the integration 
of its right-hand side in t. 

Summarising, we have cast the straight-gradient dimer (21) with > 0 and < 1 in the form of a Hamiltonian 
system 


P = 


dH 


P 

5p ’ 


9 = 


dH 


Pg = JJL 

® 89 ’ 


( 53 ) 


with the Hamilton function 


i7(p, 0, Pg) = -270 - — [Ar cosh(APe) - ax]. 


Here x = x{p,9) and r = r{p,0) are as in (47l and (49), respectively; the variable Y in (49) is defined in (45). 
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Assuming that > 1 and defining 


B. Elliptic case: |(t| > 1 


- 1 > 0 , 


the identity = 7Z^ becomes an ordinary differential equation for TZ{V): 


a \2 fdn 




(54) 


Here we used (43), while the variable r was introduced differently from (49): 


=- 


The separable equation (54) is solved by taking 


7^ = -^cos(H'P) + 


(55) 


hence 


Z = r sin(f77^). 


(56) 


This provides a simple expression for the canonical momentum: Pg =V = ^ arcsin(Z/r). 

Thus, the straight-gradient dimer (21) with > 0 and cr^ > 1 is cast in the form (1^ with 


H{p, 6, Pg) = -2"f9 -I- — [nrcos(flPe) - crx], 


where x = x{p, 9) and r = r(p, 9). 


C. Parabolic case: IctI = 1 


Assuming a — ±1, the identity +Y^ + Z^ = 7l^ reduces to an equation 

2 


on 




(57) 


with a solution 


R = r(p^ + i + 5'- 


(58) 


The rule (43) gives then Z = axV, so that 


V = a — . 

X 


(59) 


The construction of the canonical variables in the a = ±1 sector is hereby complete. In the canonical equations 


(53), the Hamiltonian is 


H{p,9,Pg) = -2j9- 


P^ + l 


y2 


® ' X2 


with X = x{p,9) and Y = Y{p,9). 
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D. One-dimensional motion: p — 0 


Finally, we consider the situation where A = B = 0 in (30a|-(30b) and hence, p = 0. Equation (31) gives 


^ /?2+4/?4 


Y = 


47^3 


U}^ 


(60) 

(61) 


The invariant manifold defined by p = 0, consists of a single curve lying in the plane (61). The parametric equations 
for this quadratic curve, with TZ as the parameter, are given by (60), (61) and (30c). 

Using the identity X'^ + Y'^ + Z'^ = TZ"^, the conserved quantity — 27 ^^^ can be represented as a function of X, Z 
and TZ'. 


-272y2 = 2'fiX^ + Z^) - 


(62) 


Here X is expressible in terms of TZ using (60) and the coefficient — 27 ^ was introduced for the later convenience. 

We choose (62) as Hamilton’s function and TZ as the canonical coordinate. The only term in (62) that remains 
independent of 7^, is 27 ^Defining the momentum by 7^ = 27 Z, the Hamilton’s function (62) becomes 


i7=^+[7(7^), [7 = 27^ 


h + 4/^4 


2 t :>2 


- cr7^ - 2YTZ: 


(63) 


while the canonical equation TZ — dH/dV reproduces equation (27). The conjugate equation V = —dH/dTZ is then 
satished automatically, because of 77 = 0 . 

Thus, the motion along the quadratic curve p = 0 is governed by a Hamiltonian system with one degree of freedom. 


VI. HAMILTONIAN STRUCTURE OF STRAIGHT-GRADIENT DIMER: <Q 

A. p > 0: two-dimensional motion 


Turning to the dimers with uP' < 0, we continue to employ the first integral p as one of the two canonical variables. 
This time, p^ is defined by (34) and equals 4|C77|. Using the explicit solution (33), one can form linear combinations 


^ , 2 - 2 / 337 ^ 47/72/73 

X+—Y = S2pe + - 

V P 2 + 4/74 

^ e , 2 - 2 / 737 ^ , 47/72/73 

^- Y = sipe + +- 

1 / P 2 + 4p4 


(64) 


where we have introduced 


sipe® = 2(7exp > ■S 2 pe ® = 277exp ) ■ 


(65) 


Provided p > 0, either of these two equations defines a real 0 which we adopt as the second canonical variable. 
Writing equations (651 in the form 


2C f2p0-vTZ 

Si — = exp --- 

p V 27 


2D fvTZ-2pe 

S 2 — = exp --- 

p V 27 


( 66 ) 


we note that since (7, 77, and p are time-independent, the argument of the exponentials in ( 66 ) is a conserved quantity. 
Therefore 


H = -2p9 -f vTZ 


(67) 


provides us with the second integral of motion for the system ( |26[ ). We will employ 77 as its Hamilton function. 

Proceeding to the construction of the momentum Pg = V canonically conjugate to 9, we note that the constancy 
of the difference lyTZ — 2p9 together with the equation (27) yield 9 = vZ. Comparing this to the canonical equation 
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(42), we conclude that the variable V should be introduced so as to satisfy the rule (43) — as in section [v| where we 
considered the case > 0. 

When was considered positive, the rule (|43|, the decomposition (|46|), and the identity =7^^ were 


the only relations necessary to derive the representations (52), (56), and (59) for the momentum V. Also used was 
the fact that x and Y were 7^-independent. When is taken to be negative, equations (64) give 


V s^e ® - sie® 47/33 

^ h - 5 - 

that is, Y remains to be 7^-independent. As for the x, we define it by 


Sie + 526 


P 2 + 4/34 


~p 


instead of (47). This preserves the validity of the decomposition (46) — and therefore, of the representations (52), 
(56), and (|59[) for the momentum Pg. As in the case > 0, equations (52), (56) and (59) pertain to \a\ smaller. 


greater and equal to 1, respectively. 

This completes the Hamiltonian formulation of the < 0)-straight gradient dimer in the part of the phase space 
with p > 0. As in the case of the > 0)-subfamily, the canonical equations are given by (53). 


B. p = 0: one-dimensional motion 


It remains to consider the invariant manifold p = 0. The manifold is described by equations (33) with C = 0 or 


D — 0 (supplemented by (30c) for the vertical coordinate). Assume, for definiteness, that D = 0. Then C and TZ 
define a pair of curvilinear coordinates on the manifold — which is, therefore, a two-dimensional surface. We will 
show that the coordinate curve corresponding to each particular value of C, is a trajectory of a Hamiltonian system 
with one degree of freedom. 


Letting D = 0, equations (33) become 


X = C exp ( 


2 - 2/337^ 
^2 + 4/34 ’ 




( 68 ) 


The second equation in (68) implies that 


H= Y + 


47 ^ 3 ^ 


exp 


27 


is a first integral of the system. We choose TZ as the canonical coordinate and appoint H as the Hamiltonian: 
H = H{TZ,'P). Here V is the momentum canonically conjugate to TZ (still to be introduced). 


The momentum should be defined so that the canonical equation TZ = dH/d'P reproduce the equation (27). The 
two equations coincide if 


Z = ;^^exp(-—7^). 




27 


(69) 


Eliminating C between two equations in (68) we can express X as 

A = -^Y + x, 

V 


where 


X = 


47/32/33 


/32 + 4/34 


— uTZ 


(70) 


(71) 


is independent of Y. Assume, for definiteness, /32/34 < 0. (The case I32Pa > 0 can be dealt with in a similar way.) 
Substituting (69) and (70) in the identity -\-Y‘^ Z'^ = TZ'^, we obtain 


Y+——x] +r] =r^, 


4/34 


dP 


(72) 
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where 


if = 


167 ^/ 32/34 


exp ( —TZ 


r -2 = (- 


V4/34 




4 / 32/34 


7^^ 


(73) 


Since neither of x, rj, or r depends on F, (72) is a separable differential equation for Y{V). A particular solution is 
Y = rsm{r]V) — :^(4/34)“^x. This relation defines V: 


V = — arcsin 
7 


*^^ 4 / 34 "' 


Lastly, we write the Hamilton function in terms of the canonical variables: 




rsiniriiV) — 4-^x ■ 

4/34 


47/33 

,/2 




where the coefficients x(72.), 77(77.), and r(77) are as in 0 and ( [73| . The two-dimensional manifold p = 0 consists of 
trajectories of the Hamiltonian system 77 = dH/dV, V = —dH/dTZ. Individual trajectories are only different in the 
value of H. 


VII. SINGULAR STRAIGHT-GRADIENT DIMER: = 0 


Finally we discuss the class of dimers with = 0. [We remind that = ^ 2(^2 + 4/34).] The straight-gradient 
dimer with /32 = 0 is gauge-equivalent to a cross-gradient system with u and v as canonical variables (section |HI[ ) . 
Therefore it remains to consider the case /32 -I- 4/34 = 0 with (32 ^ 0 only. 

We start with uncovering the Hamiltonian structure of the singular dimer with an arbitrary coefficient (33 (subsection 
VIIA). In the special case where P 3 = 0, the dimer admits an alternative, coexisting, Hamiltonian formulation. This 


is considered in subsection IVH HI 


A. General singular dimer: arbitrary value of (Ss 


The general solution of equations (28) with w = 0 has the form 

( 32(33 


X = -^YoU + ^772 - + An 


27 


(32_ 

472 


1272 

1 ~ ^3^2 


Y = --n + ^n^ + Yo, 
7 27 


(74a) 

(74b) 


where Ag and Yq are constants of integration. We define our first canonical variable by 


27 = r +-77 - ^77^ 
7 27 


(75) 


According to (74b), y is conserved: y = Yo. As in the nonsingular situation, the advantage of using the first integral 


as a canonical coordinate is that the associated momentum Py drops out of the dynamics. 

Appointing 77 as the second canonical coordinate, we need to determine the expression for the momentum P-ji 
canonically conjugate to 77. To simplify the notation, we denote it V. Note that A is not an independent variable 
here; instead, A = X{'y,TZ,V). Also note the expression for the Y component of the Stokes vector: 


r = 27 --77 +^ 772 . 

7 27 


(76) 


We define the Hamiltonian by 




f^'u3 

1272 


(77) 
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The fact that 77 is a conserved quantity follows from equation (74a): H = Xq. To identify the momentum P-ji = V, 
we compare the Hamilton equation 


7^ = 


dH 

W 


to equation (27). Since dH/&P = dXjdV, this comparison yields 

27 dV 


(78) 


Substituting (76) and (78) in the identity X'^ + = 77.^, we obtain an ordinary differential equation for X{V): 


X^ 


1 /dX 


4^2 Q'p 


(79) 


where we have introduced 


r{y,'Jl) 



y--n+^TV 

7 27 


2 


Note that r is independent of V. 
A particular solution of (79) is 


X = rsin{2'yV)\ 


then 


Z = rcos[2'yV). 


These equations define the momentum P-ji: 


P-ji =V = — arctan 
27 



(80) 


On the other hand, the momentum Py is defined by the canonical equation 


P, 


dH 
dy ■ 


Since dH/dPy = 0, the right-hand side is independent of Py and the momentum is recovered by a simple integration: 
Py = -J{dH/dy)dt. 

To complete the identification of the Hamiltonian structure of the singular dimer, we express the Hamilton function 
in canonical variables: 


132. 


27 


133. 


H{y,TZ,V) = {y-—n+^n^ + sm{2-fP) \ TZ^ - 


27 67 


y — Ti- 

7 


^^2 

27 


B. Special singular dimer: da = 0 

In this subsection, we cons ider a special subclass of sin gula r dimers where ds = 0 is satisfied along with /32 + 4/34 = 0. 
Setting ds = 0 and using (74b) to eliminate Ig from (74a), the solution (74) becomes 


02 9 02 

X = -07^^-gF7^ + Xo, 


Y = -^n + Yo. 

7 


(81) 


This time, we choose Hamilton’s function 77 to be a multiple of Ybj 

77 = -7^ - 7^, 


(82) 
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and define a canonical coordinate p as a quadratic combination of X and Y: 


p = (32Y^ - 4X. 


(83) 


Using (81), one readily verifies that p is a first integral: p = /32Yq — 4Xo. 

Appointing Y as the second canonical coordinate, the associated momentum Py = V should be introduced so 
as to satisfy the canonical equation Y = dH/dP. The conjugate equation V = —dH/dY will then be satisfied 
automatically. Comparing Y = dH/dV to and making use of (|^ gives 


dn 

W 


= 2Z. 


With the help of equation (84), the identity P? — = X^ + Y^ becomes 

1 fdP^ ^ 


^ 4 i au * ^ ’ 


(84) 


(85) 


where r stands for \/A^A|Up 2 . Using (83), X can be expressed in terms of the independent coordinates Y and p. 
This means that r is a function of Y and p — but does not depend on V : 

r' = ^(/32U2-p)2 + rU 


Accordingly, equation (85) can be considered as a differential equation for P{V). 
A simple solution to this separable equation is 


Equation (84) gives then 


P = r cosh( 2 U). 

( 86 a) 

Z = r sinh(2U). 

( 86 b) 

arctanh {Z(P). 



The relations ( 86 ) define the momentum: V = 4arctanh (Z/7^). 

Thus the singular dimer with /?3 = 0 and P 2 + 4/34 = 0 is a Hamiltonian system with the Hamilton function 

H = -yU - cosh(2U) _ p)2 + 15^2^ 


canonical coordinates p and Y, canonical momentum Py = P defined by ( 86 ), and the momentum Pp recoverable 
from Pp = —dH/dp. This Hamiltonian formulation coexists with the formulation derived in the previous subsection. 


VIII. NONLINEARITY-INDUCED PT-SYMMETRY RESTORATION 


That some conservative systems have all their trajectories confined to a finite part of the phase space, is a common 
knowledge. The harmonic oscillator provides a textbook example of this behaviour. 

Systems with balanced gain and loss may have a similar property. Assume, for instance, that the amplitudes u and 
u in (13) and (21) are small. Then these UT-symmetric dimers reduce to a two-site linear Schrodinger equation: 


iiiY V = i'yu, iv + u = —i'yv. 


When 7 is small, all solutions to this system are bounded but as 7 exceeds the critical value of 7c = 1, generic initial 
conditions lead to solutions that grow exponentially (with the growth rate A = -^ 7 ^ — 1 >0). It is common to say 
that the PP symmetry is spontaneously broken in the domain 7 > 7c (where solutions blow up) and unbroken in 
the region 7 < 7c (where all trajectories are confined). The system is said to undergo the UT-symmetry breaking 
transition as 7 is raised through 7c [IIHS]. 

Adding nonlinear terms may bring about a variety of effects. Thus, the on-site nonlinearity of the standard dimer 
([^ promotes the blow-up. In this system, large enough initial conditions trigger exponential growth regardless of the 
value of 7 [331134j : furthermore, when 7 > 1, all generic initial conditions blow up [341136] . In contrast, the nonlinear 
coupling of the cross-gradient dimer (§ softens the symmetry-breaking transition. In this case stable bounded 
solutions persist for arbitrarily large values of the gain-loss coefficient |T7]. (A similar effect is exhibited by couples 
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of damped-antidamped anharmonic oscillators |38j and solitons in a defocusing nonlinear trap with symmetrically 
distributed gain and loss [39]; hence the nonlinear softening is a general phenomenon not limited to dimers.) 

In what follows, we show that there are several classes of cross-gradient and straight-gradient T^T-symmetric 
dimers that confine all their trajectories — regardless of the value of the gain-loss parameter 7 . In these cases the 
nonlinearity not just softens the symmetry-breaking transition but suppresses it completely. The T^T-symmetry 
becomes spontaneously restored. 

The spontaneous symmetry restoration employs the same mechanism as the transition softening — just in a more 
efficient way. The exponential growth of small initial conditions is curbed by the nonlinear coupling which diverts 
increasingly large amounts of energy from the gaining to the losing site. As a result, the blow-up is arrested and all 
trajectories remain trapped in a finite part of the phase space. (Previously, a similar blow-up suppression was observed 
in a damped-driven dimer without the VT symmetry, namely, in the actively coupled waveguide pair [MlHDj .j 


A. Cross-gradient dimer 


We start with the family of the cross-gradient dimers Q. Each member of the family has two independent constants 
of motion, X and ij, where the Hamiltonian (12) has the following expression in terms of the Stokes variables: 


Jo = aiTZ^ + (x - yT. 


(87) 


Assume, first, that 0:2 < 0 while ai -b ^02 0. Noting that |y| < TZ, we obtain a lower bound for ij: 


> ( ai - ^ 


7^ 


asX — 7 — I 
2 (q:i -I- 70:2). 


(“4 + x) " 


{asX - 7 - 1 )^ 
a2 -b 4ai 


( 88 ) 


If ai is greater than 7 |q; 2 |, then, keeping in mind that X is a first integral, this inequality implies that TZ is bounded 
from above by two constants of motion. That is, there exists an TZq such that TZ{t) < TZq for all t > 0 — the trajectory 
is trapped in a finite part of the phase space. 

To find the Hamiltonian’s lower bound in the situation where 0:2 > 0, we first write Eq.(|87|) in the form 


Jo = ai 77. -b 


asX-l 

2ai 




^2 _ (Q3-A - 1)^ _ 7 

4a 1 


2 

a 2 


This representation is valid if neither ai nor a 2 is zero. Assuming 02 > 0, this gives a lower bound different from 

pj): 


Jo > ai 77 -b 


aaX-l 
2 ai 


-b 04^4 - 


2 (asA - 1)2 72 


Q !2 


4ai 


If, in addition, ai > 0, this inequality implies that there is 77o such that 77(t) < TZq for all t > 0. 

If a 2 < 0 and ai < 0, or if 02 > 0 and ai < — 7 a 2 , we can establish the boundedness of 77(t) by considering the 
lower bound for the integral —f) instead of Jo. 

In summary, trajectories of the cross-gradient dimer are confined if either (a) ai and a 2 are both nonzero and have 
the same sign; (b) ai and 02 are both nonzero and of the opposite sign, with |a 2 | < 4|ai|; (c) ai ^ 0 while a 2 = 0. 


B. Straight-gradient dimer 


Similar analysis can be carried out for the straight-gradient dimer (21). Let, first, > 0 and assume that 77 
as t tends to infinity or approaches some finite value to- Equation (30a) gives 


00 


X = —aTZ -b 0(77°) as 77 — 00 . 

This is only consistent with the inequality |X| < 77 if |cr| < 1. Consequently, if \a\ > 1, all trajectories of the 
straight-gradient dimer have to be confined: 77(7) < 77o with some hnite TZq. 

This conclusion is illustrated by the left column of Figwhich shows the surface (31) with |cr| < 1, |cr| > 1, and 
\a\ = 1. The surface is only seen to be compact in the middle panel, where \a\ > 1. 
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FIG. 2. (a) The surface (89 1 with / 32/33 7 ^ 0. In this plot, 7 = 1, /32 = 4, pz 

trajectories of the system (26l. (b) The surface (74a I with /Is = 0. Here, 7 
carries only one trajectory of the system (26 1 (depicted). 


1 , P 4 = — 1 , and Xo = 3. Also shown are two 
1, /?2 = 4, Xo = 3 and Vb = 2. This surface 


Turning to the < 0 subfamily and assuming 7?. —)■ 00 , Eq.( 33a| ) indicates that the coordinate X will grow 
exponentially in 7^ if C is nonzero. This is clearly inconsistent with \X\ < TZ. The only trajectories in ( |33a[ ) that are 
consistent with this inequality, are those with (7 = 0; here the growth becomes linear in 7?. as 7?. —>■ 00 . However if 
|cr| > 1, the slope of the asymptote of X = X{TZ) will be greater than 1. This is, again, incompatible with |X| < TZ. 
Therefore in the case \a\ > 1 all trajectories have to be bounded: TZ{t) < TZq. 

The surfaces (311 with < 0 are plotted in Fig[H right column. As in the left column, here we illustrated \a\ < 1, 

(mi 


\a\ > 1, and |cr| = 1. Only the surface with |(t| > 1 (middle panel) is compact. 

It remains to consider the situation = 0 which consists of two cases: (a) /32 = 0 and (b) P 2 + 4/34 = 0. In case 
(a), equations (28) yield 


X = Xn 


Y ^ (fe + 4/34 )Xo - 2 (33.^2 

27 27 


As 7?. —>■ 00 , the expression for Y is only consistent with |y| < TZ if Ps = 0. Therefore if ^3 ^ 0, the trajectories have 
to be confined. This is in agreement with the fact that the straight-gradient dimer with ,82 = 0 is gauge-equivalent to 
a cross-gradient dimer with ai = \p 3 and 02 = 0 , while the cross-gradient dimer with 0:2 = 0 and oi 7 ^ 0 has been 
shown to suppress the blowup (see the previous subsection). 

Finally, we let P 2 + 4,04 = 0. To avoid the duplication of results of the previous paragraph, we also require 
/32 7 ^ 0. Assuming 7^ —>■ 00 , the exact solution (74) indicates that the X component would have to grow cubically 
or quadratically in TZ if P 3 ^ Q or /33 = 0, respectively. However, neither cubic nor quadratic growth is consistent 
with the inequality |A| < TZ. Hence in the case where P 2 + 4/34 = 0 but (32 ^ 0, all trajectories have to be bounded: 
TZit) < TZq. 

To illu strat e this conc lusion geometrically, we depict the solution surface for the case (32 + 4:(3a = 0. Eliminating Yq 
between (74a) and (74b) gives 


^ = ^0 

67^ 47^ 27 


(89) 


For each Xq , equation (89) describes a surface in the (A, T, Z) space which hosts a one-parameter family of trajectories. 
If /33 7 ^ 0 and P 2 ^ 0, the surface is compact [Fig[^(a)] so all trajectories are confined to a finite part of the space. 


When (33 = 0, the surface (89) with a sufficiently large 7 (7 > |) is noncompact; yet the self-trapping of trajectories 


can be illustrated in this case as well. To this end, we note that the equation (74a) also defines a surface in the 
{X,Y,Z) space. For the given Xq and Yq, this surface hosts a single trajectory of the dimer. When (33 = 0 (while 
(^2 7 ^ 0), the surface (74aI is compact — for any 7 and any choice of Xq and Yq. This is illustrated in Fig[^b). 

To summarise, trajectories of the straight-gradient dimer are confined if either (a) ( 32{(32 + 4/33) 7^ 0 and 2|/33| > 
1/32 + 4 / 341 ; (b) (32 = 0 but (33 ^ 0 ; or (c) (32 + 4/34 = 0 but (32 7 ^ 0. 
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IX. A NOTE ON APPLICATIONS 


The standard dimer Q describes the coupling of a Kerr optical waveguide with loss to a twin waveguide characterised 
by the optical gain of equal rate [HElIloHls]. In the boson condensation context, the same pair of equations governs 
condensates in two identical potential wells [1], with one well losing and the other one being fed with atoms mm]. 
The cross-gradient dimer Q models the same potential-well geometry of the condensate, the difference from the 
standard dimer being that this time u and v are the amplitudes of the symmetric and antisymmetric state [5] rather 
than the amplitudes of the condensate in the left and right well. 

The aim of this section is to emphasise that other nonlinear models in (13) and (21) are not physically irrelevant 
either. In particular, these dimers furnish amplitude equations for couples of oscillators with physically realistic 
nonlinearities. We exemplify this correspondence by simple systems with the gain-loss balance of.two different types. 


A. Two pendula with periodic coupling 


The first system consists of two pendula with a periodically varied coupling: 

(Px . , \ Py ■ / N ^ 

-I- sin a: + pTjy = 0, -|- smw -|- k(t]x = 0. 


(90) 


The coupling is assumed to be weak and varied at the frequency close to the double natural frequency of each 
pendulum: 

K = 2 e^ cos( 2 wr), w = 1 — VtP. 

Here is a small parameter that sets the scale of the amplitude of the coupling modulation, while the coefficient 
n = 0 ( 1 ) measures the detuning of the driving half-frequency from the frequency of the linear oscillations. 

This type of parametric driving can be easily realised experimentally. For example, the pendula can be hung from 
a common horizontal rope, with a periodically varied rope tension. 

Assuming that the pendula are performing small-amplitude librations, we expand x and y in powers of e: 

x = exi + e^X 3 + ..., y = eyi + e^ys + ... 

The coefficients of the expansion are allowed to depend on a hierarchy of time scales T„ = e^oJT, n = 0,2,.... The 
times T„ become independent as e —>■ 0; in this limit, d/dr = (jj{Dq + €^D 2 -!-■■■), where = 9/9T„. We also develop 
since and siny in powers of their arguments. 

Substituting these expansions in (90) and equating coefficients of like powers of e, the order yields 

aq = -f yi = 

where A and B are functions of r 2 , T 4 ,... — but not of Tq. At the order P we obtain a pair of equations 

(£>q -I- l)x3 = 2^Dlxi - 2 D 0 D 2 X 1 - 2cos{2To)yi + ^xl, 

{Dq + l)y 3 = 2nDlyi - 2DoD2yi - 2cos(2To)xi -f ^yf. 


6 ‘ 


Substituting for xi and yi, and setting the secular term equal to zero we arrive at 

2iD2A + 2HA - iB* - ^\A\^A = 0, 
2iD2B + 2nB - iA* - ^\B\^B = 0. 
Assuming, for definiteness, > 0 and letting 

A = 2Vn{u + v), B = 2Vniu* -V*), 


(91) 


equations (91) become 


iu + v — iju = Pv* + (2|uf -I- |nf )n, 
ii) + u + pv = v^u* + (2|r;|^ -I- \u\^)u, 


(92) 


where we have introduced 7 = (2H) ^ and the overdot indicates differentiation with respect io t = Q,T 2 . The system 
(92) is nothing but the cross-gradient T^T-symmetric dimer (13) with 02 = eta = 0 and ai = 0:4 = ^. 
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B. Damped-antidamped oscillator couples 


Another, unrelated, interpretation of equations ( |13[ ) and (21) is that of amplitude equations for a damped oscillator 
coupled to an oscillator with negative damping. (This time the coupling coefficient is assumed to be constant.) 
Following the multiple-scale procedure of the previous subsection, one can readily check that the cross-gradient dimer 
(13) governs the oscillation amplitudes of the damped-antidamped pair 


Xrr + + X -\- ny = co{x^ + 5y^)x + (cix^ -I- C 2 y^)y, 

yrr -r]yT + y + KX = co{y^ + 3x^)y + {ciy^ + C2x'^)x. 


(93) 


Here k = and rj = are the coupling and the gain-loss coefficient (assumed small), while the nonlinearity 
coefficients cq, ci and C 2 can be chosen arbitrarily. The dimer parameters in (13) are expressible through the coefficients 
in (93): ai = |c2, 0:2 = ci — 3 c 2 , 0:3 = 3co, and a 4 = |ci. 

On the other hand, the straight-gradient dimer serves as the amplitude system for the oscillator couple 


Xrr + yxr + X -\- Ky = [cix^ + C2y^)x + co{ 3 x^ -l- y^)y, 
yrr -r]yr + y + KX= (ciy^ -I- C2X^)y + Co(3?/^ -I- x‘^)x. 


This time, the relation between the coefficients of the oscillators and the dimer parameters is as follows: /3i = |ci, 
P 2 = C 2 - 3ci, Ps = 3co, and P 4 = ^ 02 . 

The damped-antidamped oscillator model has been employed to interpret experiments in systems as diverse as 
a tied pair of magnetically kicked pendula [4T], two connected optical whispering galleries nails], or a tandem of 
inductively coupled active LRC circuits — one with amplification and the other one attenuated at an equal rate [43) . 


X. CONCLUSIONS 


The main results of this study can be summarised as follows. 

(1) We have introduced a four-parameter (ai, a 2 , 0 : 3 , 0 ^ 4 ) family of Hamiltonian T^T-symmetric dimers with a cross¬ 
gradient canonical structure, equations (13). The entire family was shown to have an additional integral of motion, 
independent of the Hamiltonian; hence each member of the family is Liouville integrable. All trajectories in these 
systems were described analytically. 

(2) We have considered a T^T-symmetric extension of a four-parameter {Pi, P 2 , Ps, Pi) family of conservative dimers 
with the straight-gradient Hamiltonian structure, equations (21). Unlike for the cross-gradient dimers, the original 
complex u,v variables do not constitute canonical coordinates for the straight-gradient family ( 21 ) — except when 
the straight-gradient dimer has a coexisting cross-gradient formulation. The three-parameter {Pi, P^, P 4 ) subfamily 
of straight-gradient dimers with P 2 = ^ was shown to admit such an alternative cross-gradient representation. 

We have identified canonical coordinates and momenta and uncovered the Hamiltonian structure for the entire four- 
parameter family of the straight-gradient VT -symmetric dimers. By establishing that each member of the family has 
an additional integral of motion, we have demonstrated that the entire family is Liouville integrable. All trajectories 
of the straight-gradient T^T-symmetric dimers were described analytically. 

(3) We have proved that the cross-gradient dimer with parameters 0 : 1^2 satisfying (a) aia 2 > 0; or (b) aia 2 < 0 
with \a 2 \ < 4|q;i|; or (c) ai 7 ^ 0 while 02 = 0 — has all trajectories bounded, irrespectively of the values of other 
parameters 03^4 or the gain-loss coefficient 7 . 

(4) We have established that regardless of the value of Pi and the gain-loss coefficient 7 , the straight-gradient dimer 
with parameters satisfying (a) P 2 {P 2 + 4 /? 3 ) ^ 0 , 2 \P'i\ > \P 2 + 4 / 34 |; or (b) P 2 = 0 , P 3 0 , with no constraints on (^ 4 ; 
or (c) P 2 + 4/34 = Q, P 2 ^ 0 , with no constraints on P 3 — has all trajectories bounded. 

(5) We have demonstrated that the amplitudes of libration of two coupled oscillators with on-site nonlinearities, 
driven by the periodic variation of their coupling coefficient, satisfy a T^T-symmetric cross-gradient dimer system. □ 


Thus, the principal mathematical conclusion is that the T^T-symmetric extensions of all conservative nonlinear 
Schrodinger dimers remain completely integrable Hamiltonian systems. On the other hand, the principal physical 
upshot is that there are broad classes of T^T-symmetric dimers that confine all their trajectories regardless of the 
value of the gain-loss parameter 7 . The T^T-symmetry, which is broken at the level of the underlying linear equation, 
becomes spontaneously restored thanks to the nonlinear coupling. 

The spontaneous 'PT-symmetry restoration may find applications in integrated optics where PT-symmetric non¬ 
linear Schrodinger dimers describe directional waveguide couplers. A nonlinear coupler composed of one core with 
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a certain amount of optical gain and another one with an equal amount of loss switches the entire power to one 
waveguide [3]. In the standard VT dimer, this power switching is accompanied by an unbounded power growth in 
one of the arms of the device — the growth not saturable by nonlinearity [I1II1I131137]. In contrast, no input can 
trigger an uncontrollable growth of optical modes in a dimer with the nonlinearly-restored VT symmetry. As a result, 
the VT symmetry restoration may represent a technological advantage. 
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